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Abstract 

Missions incorporated into NASA’s Vision for Space Exploration include many different destinations and 
regions; are challenging to plan; and need new and innovative trajectory design methods to enable them. By 
combining proven methods with chaos dynamics, exploration goals that require maximum payload mass or 
minimum duration can be achieved. The implementation of these innovative methods, such as weak 
stability boundaries, has altered NASA’s approach to meet exploration challenges and is described to show 
how exploration goals may be met in the next decade. 

With knowledge that various perturbations play a significant role, the mission designer must rely on both 
traditional design strategies as well as these innovative methods. Over the past decades, improvements have 
been made that would at first glance seem dramatic. This paper provides a brief narrative on how a 
fundamental shift has occurred and how chaos dynamics improve the design of exploration missions with 
complex constraints. 
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Agenda 


'/ Vision for Space Exploration 

• Exploration Overview 

* Goals and Missions 




V Traditional and Innovative Design Methods 

• A Historical Perspective of Supported Missions 

• Implementation of Chaos 

V Mission Applications 

, ■ • Sun-Earth Libration Missions 

• Earth-moon Co-linear Libration Missions 
. • Lunar Mission Design 


Vision for Space Exploration, A "Current" View 


Exploration Systems Mission 
Directorate (ESMD) 

develop a constellation of new capabilities, 
supporting technologies, and foundational research 
that enables sustained and affordable human and 
robotic exploration.” 

Themes: 

• Constellation Systems 
: * Crew Exploration Vehicle (CEV) 

- Development and Launch Vehicles 
fav: '.Exploration Systems Research and Technology 
.]■- • Prometheus Nuclear Systems 
- • Technology and Human Systems 

vAlsso jpart of the ESMD is the Robotic Lunar 
'Exploration Program (RLEP) 

Vi? L .* Lunar Reconnaissance' Orbiter (LRO) 
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Developing the vehicles and infra- 
structure that will allow us to travel 
to and explore the solar system. 



Developing the technologies today 
for tomorrow's exploration of the . 
solar system. » . 



Vision for Space Exploration, A 


‘Current" View 


Science Mission Directorate (SMD) 

Combines former enterprises of the Space Sciences and Earth Sciences 

V Solar System Exploration (SSE) (includes the former Moon and Mars 
exploration) 

X Earth-Sun System (Sun-Earth-Connections and Earth Sciences) 

X Universe (includes Origins and Structure & Evolution of Universe) 



Examples 

• Space interferometry missions, 

W-: • James Webb Space Telescope (JWST) 

* Micro Arcsecond X-ray Imaging Mission (MAXIM) Concept 

~ • SIRA(?) MMS(?) 

Space Operations Mission Directorate (SOMD) 

C X Shuttle and ISS activities 

W0T- ^ Space communications systems and the supporting infrastructure. 

X Ongoing libration orbit missions such as SOHO and WIND missions 
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An Innovative Design 


Definition: This new theory is defined as the use of 
chaos to design trajectories and orbits that can be 
used to meet complex mission goals 

Benefits: 

o Minimizes fuel cost (related to Delta-V cost) 
o Optimizes trajectory profiles 
o Provides non-standard and new orbit designs 
o Mitigates operational risks 

Other ‘synonymous’ terms 
„ o Dynamical Systems 

. -i - o Invariant Manifolds 

fv : o Capture Orbits 

yf.. o Ballistic Orbits 




A Sample, of Analysis 

So let us look at a few sample ESMD 
and SMD missions: 


Sun-Earth libration orbits 
Earth-moon libration orbits 
Lunar mission design 
Use of chaos to aid in their design 


Historical and Future Missions in Libration Orbits 



COMET HAUW 

>28- es 


HALO ORBIT 
6 uarBAva 
5 WORST 


ISEE3 MANEUVERS FROM LAUNCH 
TO HALO ORBIT 
TO COMET EXPLORATION 


Mission 


Location/ 

Type 


Amplitudes 

j(Ax,Ay,Az) 


Launch ^Total AV- 
' Year : Allocation 

1978 430 


Transfer-, 


LlHalo/L2/Comet 175000, 660670, 120000 

l* 1 mission ■ ' 1 ■ 

Ll-Lissajous 10000, 350000, 250000 


Direct 


Multiple 
Lunar Gravity 
Assist 
Direct 


WIND' 


SOHO 


666672, 120000 


Ll — Lissajous 
1** small amplitude 
L2-Lissajous 
1* L2 Mission 
Ll-Lissajous 


Direct . 
(Constrained) 
Single Lunar 
Gravity Assist 
Direct 


81775, -264071, 157406 


n/a, 2 64000, 264000 


250000, 800000, 250000 


Genesis 


Ll-Lissajous 
Launch 
Constrained 
L2-Quasi-Periodic 
Lissajous 
L2-Lissajous 
Tethered 
Formation 
L2 - Lissajous 
Formation 
L2 - Lissajous 
Loose Formation 
Ll-Lissajous 
Large Lissajous 
L2 - Lissajous 
~30 S/C Formation 

L2 - Lissajous 
Formation? 


JWST 


290000, 800000, 131000 


Direct 


Direct 


290000, 800000, 131000 


Large Lissajous 


Direct 


MAXIM 


Large Lissjaous 


Single Lunar 
Gravity Assist 


Constellatioa-X 


300000, 800000, 350000 


Darwin 


Direct 


Large Lissajous 


Stellar Imager 


Lissajous 



ISEE-3 / ICE 


t“ Halo Orbit 


If 


I . HALO-ORBIT INSERTION 

\ i 

\ l 

\ r" • HALO ORBIT 

' | PER 100 v I MONTHS) 


LAUNCH C, '--0 Jt Km>/ W 


_,unar Orbit 


^Mission: 


‘ Ljssajous Orbit: 




Solar-Rotating Coordinates 
Ecliptic Plane Projection 

Investigate Solar-Terrestrial relationships, Solar Wind, Magnetosphere, 
and Cosmic Rays 

September, 1978, Comet Encounter Sept., 1985 

LI Libra tion Halo Orbit, Ax=~l 75,000km, Ay = 660,000km, Az~ 

120,000km, Class I 

Mass=480Kg, Spin stabilized. 

First Ever Libration Orbiter, First Ever Comet Encounter 

Fanjntar ct al [I9S5] Trajectories and Orbital Maneuvers Tor the 1SEE-J-1CE. Comet Mission. JAS 33. No. 3 


;i||ajoiis Orbit 
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GEOMAGNETIC TAn. 


Solar-Rotating Coordinates 
Ecliptic Plane Projection 


Investigate Solar-Terrestrial Relationships, Solar Wind, Magnetosphere 

November, 1994, Multiple Lunar Gravity Assist 

Originally an Ll Lissajous Constrained Orbit, Ax~10, 000km, Ay ~ 

350,000km, Az~ 250,000km, Class I 

Mass=1254kg, Spin Stabilized, 

First Ever Multiple Gravity Assist Towards Ll 









Mission: ”4 

hch: >" 
Lissajous Orbit: 


cef raft: 




Produce an Accurate Full-sky Map of the Cosmic Microwave Background 
Temperature Fluctuations (Anisotropy) 

Summer 2001, Gravity Assist Transfer 

L2 Lissajous. Constrained Orbit Ay~ 264,000km, Ax~tbd, Ay~ 264,000km, 
Class II 

Mass=818kg, Three Axis Stabilized, 

First Gravity Assisted Constrained L2 Lissajous Orbit 


JWST 


1 l- fiKlV : :ia : . t f/C/ II 

471.28SB days 


jit- f iTes.. - ■ Uicus-v <Last> 


E3rt h 'I \ jj 


v|^M|ssion: 


■fast 


aunch 

~.«S- - . 

bus Orbit: 


Solar-Rotating Coordinates 
Ecliptic Plane Projection 

JWST Is Part of Origins Program. Designed to Be the Successor to the 
Hubble Space Telescope. JWST Observations in the Infrared Part of the 
Spectrum. 

JWST-2012, Direct Transfer 

L2 large lissajous, Ay- 294,000km, Ax-800, 000km, Az~ 131,000km, Class I 
or II 

Mass~6000kg, Three Axis Stabilized, ‘Star’ Pointing 

Observations in the Infrared Part of the Spectrum. Important That the 

Telescope Be Kept at Low T emperatures, ~3°K. Large Solar Shade/Solar Sail 
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A General Design Process 

In Low Fidelity Software 

• Use Chaos mathematical expressions for preliminary orbit design, e.g. Circular 
Restricted Three-Body (CRTB) problem. 

• Generate orbit families via differential correction and continuation. 

• Analyze the properties of these orbits and to meet mission requirements. 

• Obtain orbit architectures. 

• Apply two-step differential correction scheme to selected orbits. 

• Add multiple revolutions for baseline mission duration. 

In Higher Fidelity Software: 

Differentially correct in full ephemeris model. 

, Constrain orbit to desired goals, apply chaos to obtain Av. 

A - ’ " Acquire Av and fuel budget for station-keeping by perturbing initial target 
% states in unstable directions and adding Av errors. 

'•.'Analyze mission requirements and constraints (e.g Sun angle limits and 
j.- Facility access) 




Chaos - System Application . 


Numerical Systems 


Chaos Systems 


IsMiiiP 


• Limited Set of Initial 

Conditions 

• Perturbation Theory 

• Single Trajectory 

• Intuitive DC Process 


• Qualitative Assessments 

• Global Solutions 

• Time Saver / Trust Results 
« Robust 

• Helps in choosing numerical 

methods 

(e.g., Hamiltonian => 

Symplectic Integration Schemes?) 





Chaos and Invariant Manifolds 

s Use of invariant manifolds are directly applicable to weak stability 
boundary and libration trajectory design 

S Together with differential corrections, the use of invariant manifolds 
provides an efficient method to obtain transfers and control 

■/ Invariant manifolds results can be used as a initial conditions for NASA 
mission design software 


Chaos System Approach Transfer 

Design a Large Libration Orbit’s Transfer Trajectory - 
Projections of All Invariant Manifolds for Time Interval 


Portion of a Manifold - PI P2rot Frame 



2 +x-axis 
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L 2 Butterfly Orbit Family 


Applications of Chaos in the Earth-Moon Region 


o A natural transfer between co-linear orbits is defined by the Eigenstracture of the co- 
linear STM and the dynamics in question. 

o The co-linear unstable and stable modes are in the familiar directions and indicate that 
direct transfers between libration points are straightforward 

o Use the unstable mode of Lj/I^ to perform a AV {possibly near or at zero magnitude) to 
transfer to the stable mode of Lj/L) once the trajectory has moved to the other side of the 
secondary mass. 

Along - | ■ - j | . , -..at ' ( 1 1 Earth-Moon Co -1 inear 

Unstable Modes Departures from L, and Lj 

Returns to ; | ** ’ • *4^ I 

i I | .„ — 

; |WPopans 1 -Afa ’ 

ij»„ . ***** i.jioj J 

illibrattoh.Pbinis ’itUt ~USSS ‘sUSS » -Uiii 
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Applications of Chaos in the Earth-Moon Region 


/ 

Manifold Tube 


Manifold Interse 


Moon 


: f lacing a facility at the South Pole of the Moon poses questions concerning the orbital architecture of 
f Vie communicating satellites. Constant communication can easily be achieved with Earth-Moon 
frtibfaiion point orbits. We analyze different architectures for nearly rectilinear halo orbits, vertical 
..orbits, and other three-body variations for lunar coverage of the South Pole. Using invariant 




Applications of Chaos in the Earth-Moon Region 


Invariant Manifolds 



• Dynamical system provides the structure 

• Numerical DC used for targeting process 

• Need to improve methods 

• Example: Finding intersections of Sun-Earth and 
Earth-Moon Manifolds for Transfer Trajectories 





Comparison of Direct and Weak Stability Boundary 
Transfer to a 100km Circular Lunar Orbit 


Direct Transfer Weak Stability Transfer 


Launch Vehicle Injection (km/s) 

3.13 

3.17 

Launch C3 (km 2 /s 2 ) 

-2.11 

-0.59 

Equivalent launch mass for C3 

. 1595 

1545 

Total Delta-V to attain mission orbit (m/s) 

821 

664 

Trip Time (days) 

.4.5 

98.0 

Max distance from Earth (million km) 

0.367 

1.5 


>AV Improvement 
of 19.12% 


,^^^^iar:bm^j|u.$sion orbit 
jriutltude . 
fh'Vdiicle of Delta-II 




m^sjoomputed via rocket equation 
|||^ngdaurich mass, lsp=220, 
above 


New Orbits and Future Challenges 

Upcoming missions also bring new challenges that individually may easily be 
met, but in combination they become problematic. These may include: 


ISP 
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Lunar Orbits for Relay Spacecraft 

Biased Orbits when using large sun shades 

Frequent Spacecraft Perturbations (momentum unloads) 

Constrained communications 

Shadow restrictions 

Very small libration orbit amplitudes (clOOOOkm) 

Limited thruster directions 

Transfers Between Libration Orbits and the Moon 

Earth-Moon libration orbits 

Continuous control to reference trajectories 

Quasi-stationary orbits 

Human exploration 

Servicing of resources in libration orbits 


WEBas** 











History, Definitions, & Modeling 


Mathematical Histor 


Poincare 


> Defined three body problem in work on ^ Stability of motion and use of potential 


lunar motion. 

> Proved existence of co-linear points 


> Development of equilibrium points 


functions 

> First to recognize the need for a qualitative 
approach to three body problem which is 
unsolvable in closed form 

Jacobi 

> One exact integral of three body system 


Definitions & Modeling 

^'Easiest to model the system as the Circular Restricted Three Body Problem 
(CRTBP) where m 1 >>m 2 >>m 3 

.^nij - primary, m 2 - secondary, m 3 - body of interest 
Emotion of Earth about Sun is circular 
H Emotion of m 3 is in plane of nr^ & m 2 
>CRTBP can be solved exactly 

^Unfortunately, unmodeled forces (solar radiation pressure, other 
I r gravitational bodies - Jupiter, etc.) and physical reality (non-circular motion 
gror EM system about sun) cause perturbations 33 


Libration Points 


- Lj and_L 2 distance of 1.5 million km 
and L s distance of 150. million km 


What Are They?? Where Are They? 

> Equilibrium or libration points represent > Collinear Points: L 1; L 2 , L 3 (unstable) 
singularities in the equations of motion > Triangular Points: L„ L s (stable) 
where velocity and acceleration 

components are zero and the forces are p 

balanced / \ 

> Viewed in the rotating frame: centrifugal / L \ 

(Coriolis-Type) force balances with / \ 

gravitational forces of the two primaries - /. 

£C / 3.U. y 

Libration points are in plane with no Z > \ 

• component. Orbits are mapped to a L 3 Suj^ L, \ L, 

rotating frame where there are no time ° ~ 

I - dependent forces \ 00, . iL "- ' 

. . / Eartb/Moon 

mmm r system of interest involves the Sun \ / System 

^^(mj), the Earth-Moon system (m 2 ) and \ / 

IpthtHe spacecraft m 3 \ / 


' T ' 
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Solar Rotating Coordinate^ 
Ecliptic Plane Projection 
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